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Abstract 

Theory of long-range interactions between an atom in its ground S state and a linear molecule 
in a degenerate state with a non-zero projection of the electronic orbital angular momentum is 
presented. It is shown how the long-range coefficients can be related to the first and second-order 
molecular properties. The expressions for the long-range coefficients are written in terms of all 
components of the static and dynamic multipole polarizability tensor, including the nonadiagonal 
terms connecting states with the opposite projection of the electronic orbital angular momentum. 
It is also shown that for the interactions of molecules in excited states that are connected to the 
ground state by multipolar transition moments additional terms in the long-range induction energy 
appear. All these theoretical developments are illustrated with the numerical results for systems 
of interest for the sympathetic cooling experiments: interactions of the ground state Rb(^S) atom 
with C0(3n), 0H(2n), NH(iA), and CR^n) and of the ground state Li(2S) atom with CR^U). 
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I. INTRODUCTION 



Recent developments in laser cooling and trapping techniques have opened the possibility 
of studying coUisional dynamics at ultralow temperatures. Atomic Bose-Einstein conden- 
sates [l| are of crucial importance in this respect since investigations of the collisions between 
ultracold atoms in the presence of a weak laser field leads to precision measurements of the 
atomic properties and interactions. Such collisions may also lead to the formation of ul- 
tracold molecules that can be used in high-resolution spectroscopic experiments to study 
inelastic and reactive processes at very low temperatures, interatomic interactions at very 
large distances including the relativistic and QED effects, or the thermodynamic properties 
of the quantum condensates of weakly interacting atoms 2j. 

Recently, experimental techniques based on the buffer gas cooling jsl or Stark deceleration 
J] produced cold molecules with a temperature well below 1 K. Optical techniques, based on 
;he laser cooling of atoms to ultralow temperatures and photoassociation to create molecules 
sl, reached temperatures of the order of a few /iK or lower. Spectacular achievements 
were reported only very recently with the Bose-Einstein condensation of homonuclear alkali 
molecules starting from fermionic atoms js Sl. 

A major objective for the present day experiments on cold molecules is to achieve quantum 
degeneracy for polar molecules. Two approaches to this problems are used: indirect methods, 
in which molecules are formed from pre-cooled atomic gases JqI, and direct methods, in which 
molecules are cooled from room temperature. The Stark deceleration and trapping methods 
pioneered by Meijer and collaborators 



10| are the best developed of the direct methods and 
provide exciting possibilities for progress towards quantum degeneracy. 

Beam deceleration can achieve temperatures around 10 mK. However, condensation re- 
quires sub-microkelvin temperatures. Finding a second-stage cooling method to bridge this 
gap is the biggest challenge facing the field. The most promising possibility is the so-called 
sympathetic cooling, in which cold molecules are introduced into an ultracold atomic gas 
and equilibrate with it. Sympathetic cooling has already been successfully used to achieve 
Fermi degeneracy in ^Li ll\ and Bose-Einstein condensation in "^^K 12|. However, it has not 
yet been attempted for molecular systems and there are many challenges to overcome. In 
Berlin, Meijer's group has now developed the capability to trap ultracold '^^Rb atoms for use 
in sympathetic cooling jl3|. In London, Tarbutt's group has begun to set up experiments 
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using an ultracold gas of ^Li atoms to cool molecules Q. Thus far, open shell molecules 
like C0(3n), 0H(2n), NH(iA), and CR(^U) could be decelerated, and are good candidates 



for sympathetic cooling. Other simple systems like LiH (15], ND3 [ly, [iTj, and ions, Yb 



and Ba^ 
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g|-|2l|. were also investigated, both exprimentally and theoretically. 
Very little is known about collisions between polar molecules and alkali metal atoms, and 
results of theoretical studies on them are essential to guide the experiments and later to 
interpret the results. There are two essential ingredients: good potential energy surfaces 
to describe the interactions, and good methods for carrying out low-energy collision cal- 
culations. Hutson and collaborators has pioneered the study of potential energy surfaces 
for interactions between polar molecules and alkali metal atoms [16, 



17, 



22 . However, the 



theoretical methods available for calculating the surfaces have some significant inadequacies, 
which need to be addressed before quantitative predictions will be possible. 

When dealing with collisions at ultra-low temperatures the accuracy of the potential in the 
long range is crucial. Therefore, the methods used in the calculations of the potential energy 
surfaces should be size-consistent 2^ in order to ensure a proper dissociation of the electronic 
states, and a proper long-range asymptotics of the potential should be imposed. The latter 
task is highly nontrivial when a molecule is in an open-shell degenerate state. To our 
knowledge this problem has only been addressed by Spelsberg [24 1 for the CO-I-OH system. 



and by Nielson et al. 



25l | and Bussery-Honvault et al. 26|, \2J\ for an atom interacting with 



open-shell molecule. The latter considerations were limited, however, to the Cg coefficients 
at most. Standard approaches based on the symmetry- adapted perturbation theory within 



the wave function 
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29l | or density functional formalisms 30| fail in this case. 



In this paper we report a theoretical study of the long-range interactions between an 
atom in its ground S state and a linear molecule in a degenerate state with the projection 
of the electronic orbital angular momentum A. The expressions derived in this work are 
applied to systems of interest for the sympathetic cooling experiments, i.e. to interactions 
of the ground state Rb(2S) atom with 00(^0), 0H(2n), NH(iA), and CH(2n), and the 
Li(^S) atom with CH(^n). The plan of this paper is as follows. In sec. [Ill we present the 
derivation of the long-range asymptotics for the interaction of a ground S state atom with 
linear molecule in a degenerate state A. We discuss which polarizability components of 
an open-shell species are needed to express the asymptotic interaction energy and how the 
expression for the interaction energy depends on the adopted basis (spherical or Cartesian). 
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In this section we also show that for molecules in an excited state that is connected to the 
ground state by multipolar transition moments, a new term in the long-range expansion 
appears. In sec. IIIII we present the computational approach adopted in the this paper, 
discuss our results, and compare them with the data from the supermolecule calculations. 
Finally, sec. [IV] concludes our paper. 



II. THEORY 



We consider the interaction of an atom A in the ground S state iPa{^) and a linear 
molecule B in a state ipsi-^), where A is the projection of the electronic orbital angular 
momentum of the molecule on the molecular axis. All the quantities relating to the atom 
and the molecule will be designated by subscripts A and B, respectively. Since we are 
interested in the long-range interactions, the resulting spin multiplicity of the complex do 
not play any role in our further developments, and will be omited. The electronic state of 
the molecule does not need to be its ground state. The Hamiltonian H of the complex AB 
can be written as: 

H = Ho + V, (1) 

where Ho is the sum of the Hamiltonians describing isolated monomers A and B, Hq = Ha + 
H^, and V is the intermolecular interaction operator collecting all Coulombic interactions 
between electrons and nuclei of the monomer A with the electrons and nuclei of the monomer 
B. Assuming that the electron clouds of the monomers do not overlap, V can be represented 
by the following multipole expansion 



34 



36|: 



V 



Ia,Ib=0 



Ia+Ib 

^-Ia-Ib'I ^ ( 
m=-lA-lB 



lA+h 



(2) 



where the constant C;^,;^ is given by 



47r 



21a + 2/b + 1 



1/2 



21a + 21b 
21a 



1/2 



(3) 



Yi"^{R) is the normalized spherical harmonic depending on the spherical angles R of the 
vector connecting the centers of mass of the monomers A and B in a space-fixed coordinate 
system, denotes the multipole moment operator in the space-fixed frame. We made also 
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use of the coupled product of two spherical tensors: 

Qi^®QiX = E {U,mA;lB,mB\L,M)Q';i^Q'^/. (4) 

where {IaiTtla,Ibi^b\L., M) is the Clebsch-Gordan coefficient. 

A state V^_b(A) of a linear molecule with A 7^ is doubly degenerate, and so is the state 
of the complex AB at -R = 00, which is just a product '?/'^(S)'?/'b(A). For the interaction of a 
ground state S atom with a molecule, the first-order electrostatic energy vanishes identically 
in the multipole approximation, so the degeneracy is lifted in the second-order, leading to 
the splitting of the two states ip a{^)'4''b{'^\^\) 1 into the A' and A" states of the complex AB 
at finite R. Thus, to obtain the long-range behavior of the A' and A" states we have to 
diagonalize the second-order interaction matrix: 

The elements V?\, of V'-^^ are given by the standard expressions of the polarization 



theory 



37| and can be decomposed into the induction and dispersion parts: 



(2) _ >^, (^A(S)V;i,(A)|V|V;^({A:})^B(±|A|))(^^({fc})^B(±|A|)|y|^A(S)^i.(AO) 

^ A,A' - 2^ - 
{k} {11} 

where Ws.fc = — is the excitation energy of the atom from the ground state ^Ja^^) 
to the excited state ipA{.{k}) characterized by the set of quantum numbers {fc}, and c^A.n = 
— E\ is the excitation energy from the state '?/'b(A) to the excited state ipB{{n}) of 
the molecule with the set of quantum numbers denoted by {n}. In the above equation the 
sign prime on the summation symbol means that the ground state of the atom or the states 
'?/'b(±|A|) of the molecule are excluded from the summations. 



It can easily be shown 



25| that in the case of interaction with an atom in S state the 



(2") (2') 

matrix elements ^ and V_l _^ are equal, and the same holds for off-diagonal elements 
_^ and V_^\. The eigenfunctions and the corresponding eigenvalues of the matrix can 
simply be constructed on the symmetry basis only, as the proper zero-order wave function 
of the complex AB should have a defined behaviour under the refiection in the plane a 
containing the three atoms: ctiJjab = i^AB ^' state and (JiPa^b — (~1)^ab ^" state. In 



our approach the three atoms he in the xz plane and the transformation between the space- 
fixed frame and the body-fixed frames is chosen in such a way that it leads to the coincidence 
of the axis y in all relevant frames of reference, see Fig. 1. Thus, the symmetry (A' or A") 
of the zero-order wave function t/'^ characterizes its behaviour under the refiection in the 
plane a^z- It is determined by the symmetry properties of the wave functions of the two 
constituent monomers, i.e. '?/'yi(S) and iPb{^), with respect to the refiection in the plane axz 
of their body-fixed frames, which happens to be the same. The symmetry relations for the 



monomer electronic wave functions are (A 7^ 0) 33|, |38|] : 



a..VB(A) = (-I)^^b(-A), a^zipA^^) = (-I)^-^^a(S), (7) 

where pa defines the spatial parity of the atomic wavefunction in the S state. Therefore, the 
properly adapted wave functions read: 

^SS = ^ [^a(S)^b(A) + (-1)^+^-^+^^^(S)^b(-A)] . (8) 

and the corresponding second-order energies (i.e. eigenvalues of the matrix V*-^-*) are 

^{S = <1 + (-1)^^"^^M'1a (9) 
with / = for A' and / = 1 for A" state. A standard phase convention of Condon and 



Shortley j40| was adopted to define V'b(^) ^-^id iPb{.~^)- More details on the symmetry 



properties of the electronic wave funcions of diatomic molecules can be found in Refs. 



38 



33 



391]. Here we only want to conclude that in the case of the systems under study the A' 
state corresponds to the combination of the two terms with the minus sing in Eqs. (IE]) and 
(E]) for Rb+0H(2n), Rb+CO(='n) and Li+CH(2n) , while the A' state has the plus sign for 
Rb-hNH(^A). The opposite holds for the A" states. 

The multipole expansion of Va,a' is readily obtained by inserting the multipole expansion 
2]) of the interaction operator V into Eq. ([6]) and collecting terms, as it was done in Refs. 
. Specifically, the derivation is based on the well known transformation properties of 
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the multipole operators from the space- fixed (with the index mx) to the body- fixed (with 
the index kx) frame of each monomer (X=A or B): 

E ^Sx(^^)^f"' (10) 

kx=-lx 
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where D'J^l'f^^{Rx) is the Wi gner D matrix, and on the addition theorems for the D 
functions and spherical harmonics: 



Lx Kx=-Lx Mx=-Lx 

^iLxr 



yi:ZBiR)-yrjijR) =E E 



X {lx, mx; I'x, m'x\Lx, Mx)D\'^^;^^{Rx), (11) 
{21a + 21B + 1){21a' + 21b' + 1)1 



47r(2L + 1) 

L M=-L \ ' J 

X {I A + -m; I A' + Ib', -m'\L, M){Ia + /b, 0; /^^ + 0|L, 0) 
xYt^{R). (12) 

The sets of the Euler angles Ra and Rb describe rotations of the space-fixed frame to 
the appropriate body-fixed frames. For our convinience we choose the Z axis along the 
intermolecular axis connecting the two centers of mass. Then, the angular factor Yl\R) 
reduces to: 

Yt{R) = yt{^A)=[^^) Sm,o. (13) 

For an atom in an S state the quantum numbers La, Ma and Ka are zero, so Dm^Ika^R^) ~ 
1. For an open-shell linear molecule we have Mb = 0, and the set of the three Euler angles 
is Rb = {0,6,0), so 



<f (0,^,0) = (-l)"-y/ |^^;gj P.?(cos^) (14) 

where P^^{cos6) denotes the associated Legendre polynomial. The set of the Euler angles 
Rb = (0, 6, 0) for a linear molecule is consistent with the foregoing requirement of the 
coincidence of the axis y in all relevant frames of reference, cf. Fig. 1. The other possible 
set would be Rb = {Sn /2,6,7t /2) which leads to the coincidence of the clXlS OC clS it was 
adopted in Ref. |25| . 

It is useful to express the final equations for the elements of the matrix V^^^ in terms of 
the static and dynamic multipole polarizabilities of the atom and molecule. For an S atom 
we use the standard definition: 

. ^ {MS)\Qr\M{k})){M{k})\Qf\MS)) ...x 

«mm'(^) = 2 > Ws,fc 2 2 ' (15) 

{k} s.fc 



while for an open-shell linear molecule we introduce extra superscripts A and A' to distinguish 
between diagonal and off-diagonal components: 



{Mmr\M{n})){M{n})\Qf\M^')) 



(16) 



The corresponding irreducible polarizabilities are obtained by Clebsch-Gordan coupling: 



(17) 



m,m' 



The only nonvanishing components of the irreducible polarizability for an atom in the S 
state are ag'^^, while for an open-shell linear molecule the nonzero components are ^'^'a^' ■'^ 
with is: = for A = A' and with if = 2A if A' = -A. 

We are now ready to give final expressions for the long-range coefficients expressed in 
terms of the irreducible components of the polarizabilities: 



T/(2) 



''''^&".^M)+^^^'c^^^.,(disp) 



Pi^'(cos( 



(18) 



where the constant f, ', ,, is given by: 

'{21a + 21b + 1)K21a + 21b' + 1] 



cL,K 
'IaIbI'b 



X 



(2/a)!(2/b)!(2/^)!(2/bO! 
/2L + 1\^/^ 1{L-K) 



^'UIa + Ib Ia + I'b L 




V 2/^ + 1 



Is I A + Ib I a 

^L + ^y-\lA + l'n I'b L 



(19) 



and the expressions in the round and curly brackets are the 3j and 6j coefficients, respec- 
tively. Combining all terms with the same power n = 2/^ + /5 + /^ + 2 in the above expansion, 
we will get standard long-range coefficients C^'^: 



by means of which the asymptotic expansion of Eq. simply reads: 



(20) 



(2) 
(±) 



EE 

71=6 L 



^P°(cos( 



(21) 
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The dipersion part ^'^'C*;^';^;/ (disp) is proportional to the Casimir-Polder integral over 
the atomic and molecular polarizabilities calculated at imaginary frequencies: 

1 /'OO 

-'^'Cf;,^,, (disp) = ^l^ a^'-n..) ^'^'ar^%u)du, (22) 

while the induction term ^'^'C*;^'/^;/ (ind) is the product of the static polarizabihty of the 
atom and permanent multipole moments of the open-shell molecule: 



^'^'cf,^, (ind) = a^'-)°(o) liMmiMMM)) ® {MMM)\Qi'MBm]'- (23) 



B 



J K 



(2) 

A few comments are needed here. The expression for the diagonal term j( is the same 
as for the interaction between atom and linear molecule in a spatially nondegenerate state 
(S). The additional term emerging when A 7^ is the off-diagonal _yv It has both 

induction and dispersion parts. As it was already stated, the only nonvanishing term in 

(2) 

y\ _A appears for K = 2A. The lowest order long-range coefficient, to which the off-diagonal 
term contributes to, is C^'^ with n = 4 + 2|A| for the dispersion part and n = 5 + 2|A| for 
the induction part. The source of the induction energy in comes from the fact that 

open-shell linear molecules have an additional independent component of the permanent 
multipole moments, {tlj{A)\Qf^\ilj{—A)), in contrast to the S state molecules with only one 
independent component, namely {ip{A)\Q^\ilj{A)) . It is obvious that the second independent 
component appears for / > 2|A|, therefore in the case of the 11 states there is no induction 
contribution to the off-diagonal Cg'^ coefficient. The dispersion terms in V^^l^ results from 
the presence of additional components of the polarizabihty terms for the open-shell linear 
molecules. 

In general, the number of independent diagonal terms ^'^a"^, is equal to 2/< + 1 (where 
/< is the smaller of I and I'). Each of the non- redundant components ^'^a^'^, comes from 
the transitions in the sum (fT6|) through the intermediate states ip{{n}) with the projection 
of the eletronic angular momentum ranging |A| — (2/< + 1) to |A| + (2/< + 1). Due to the 
transformation properties, the number of independent irreducible components ^'^a^^ will 
also be 2/< + 1. 

In addition, for open-shell linear molecules, there are nondiagonal components of the 
polarizabihty tensor ^''^a^'k' which do not vanish, and are not related to the diagonal 
terms. The nondiagonal terms in the multipole polarizabihty tensor appear if the condition 



k' + k = 2A can be fulfilled for a given set of the quantum numbers {l,k) and {l',k'). 
They result from the fact that there are transitions through the same operator to states S"*" 
and S~, and then ^'^a^'k the diffrence between these two separate contributions (e.g. 
^'^<yl\ for the n state). The other mechanism leading to the appearance of the nondiagonal 
^'"^'^fcfc' terms is parallel with the source of the off-diagonal induction terms, namely there 
are possible two independent transition moments between the state in question \ip{A)) and 
the exited states through an operator of the same rank, e.g. ^'~^aQ2 the 11 state, to be 
compared with the diagonal ^'^a^ . 

Let us discuss the dipole polarizabilities in more details. For open-shell diatomic 
molecules there are three independent spherical components ^'^a^^/ with {k, k') equal to 
(0,0), (1,-1) and (—1,1), the corresponding transitions in sum-over-states occur through 
the excited states with a projection of electronic angular momentum equal to A, A — 1, and 
A + 1, respectively. In case of the 11 states there is extra nondiagonal component with 
intermediate states S"*" and S~ in the sum-over-state expression, which come with opposite 
sign, in contrast to ^'^a^\i in which states S"*" and contribute with the same sign. The 
four corresponding irreducible dipole polarizability components would be ^'^Oq^^'*'^, ^'^Oq^^^^, 
^'^Oq"'^^''^, and ^'"^^2^"'^^^, the last one is nonvanishing only for the 11 states. 

From the analysis of Eq. (fT8|) it turns out that not all nonvanishing irreducible compo- 
nents ^'^'a^' ■'^ of the molecule polarizability tensor are needed to express the interaction 
energy of an open-shell linear molecule with an atom. The dipole component ^'^o^q^^^ is 
redundant in this case. It follows from the expression ( |T9|) that all terms aj^ with 
l + V + L odd do not contribute to the interaction energy in the second order E(^j^y However, 
if we had an asymetric molecule or an other linear open-shell species instead of an atom then 
all non-zero irreducible componets would be necessary. For a A state there are three inde- 
pendent dipole polarizabilities, but again, in the case of interaction with an atom, ^'^Oq^^^^ is 
redundant. The first off-diagonal term for a A state will be '^~'^ol\^^'^ or equivalently 

Sometimes it is more convenient to use the Cartesian basis both for the states |A) and 
multipole moments. The transformation formulas between the two bases can be found in Ref. 
[sol . We focus again on the 11 states. In the Cartesian basis the two degenerate states | ± 1) 
are usually referred to as jllx) and jlly). The four independent dipole polarizability tensor 
components would be ""'""axx, "^'""ci^yy) ""'""ttzz, and "^'^Oxy, where we adopted the index (x,y) 
in place of (-1,1) to distinguish between particular diagonal and off-diagonal components. 
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It is possible to define the fifth Cartesian component, namely ^'^Oyx, however, due to the 
relation: "^'""axx — "^'"^ctyy = "^'^Oixy + '''•^ttyx, it will not be independent |2J]. To see better why 
the Cartesian basis may be useful let us mention how the irreducible spherical components 
are related to the Cartesian ones for the 11 states: 

1,1 (1,1)0 -1-1 (1,1)0 1 /x,X^ I X,X^ I X,X^ N 

"O = = ( ttxx + ttyy + ttzz) 



1 
V6 



=-1.-1 = — (2 ^'^a,, - ^-^a^ - ^'^a^y, 



''-'a^''^' a^T = ^'^"yy - -cxx (24) 

An inspection of the above relations shows that in order to express the interaction energy 
of a n molecule with an atom, only diagonal Cartesian components are needed, as the 
off-diagonal ^'^a^y contributes only to ^■^aQ^'^''^, i.e. the term not present in expression for 
E^^y The off-diagonal irreducible component ^'^^0:2^'^^^ is equal to the difference between two 
perpendicular a's calculated for one of the Cartesian states. This observation is more general 
and it turns out that if we decide for a Cartesian representation of the |A) states, then all off- 
diagonal irreducible components ^'^^agl^^ can be related to some diagonal Cartesian terms. 
Still, we will have to calculate some additional Cartesian components, however, diagonal 
only, which are not present in the expression for the diagonal ^'^o'q terms. Calculations 
of the off-diagonal polarizability components in the Cartesian basis are indispensable if the 
interacting system comprises of an open-shell molecule and asymmetric species. It means 
that in such a case the interaction energy in the multipole approximation is expressed in 
terms of, for instance, ^'^cx^y component. It is worth noting that ^'^cx^y is irrelevant in the 
describtion of the interaction of the open-shell molecule with the extrenal electric field in the 
Stark effect. Thus, the induction energy and, consequently, intermolecular forces explores 
properties of the open-shell diatomic molecules which are not accessible otherwise. 

Note that Eq. fl2^ is strictly valid only when the molecule is in its ground electronic state. 
If the molecule is in an excited state that is connected to the ground state (or to any other 
state lower in energy) by multipolar transition moments then the Casimir-Polder integral is 
no longer valid and an extra term has to be added to the energy. The reason behind this is 
the property of the Casimir-Polder integral that if the two elements in denominator, and 
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eb, have opposite signs then we obtain (assuming that < 0): 

2 r eA es , 1 



9 o 9 = , — T, (25) 

^ Jo ^A + '^ <^B+'^ eyi + lesl 

instead of the value of l/{eA — \^b\), which we want to decompose into the product of two 
terms depending on monomer properties only. Formally, we may write the following identity 
(again < 0): 

1 2 r eA , , 1 



, 2 9 9 9^^ + , — r + , — r- (26) 

CA-Iesl vr Jo (^a + ^ + ^ + \eB\ - jesl 

For any molecular state with < we have the summation over all atomic states with 
positive excitation energy e^, this means that the two last factors in the above equation will 
add up to yield dynamic polarizability of atom at frequency u = I^bI- Therefore, if we want 

(2) 

to express the whole interaction energy in the second order E^^^^ in terms of the monomers 
properties only, then we are forced to add an extra term to the dispersion part, Eq. fl22l) . 
in order to compensate an error introduced by the integral representation of the dispersion 
energy. This additional term depends on the dynamic polarizabilites of the atom calculated 
at frequency u equal to the energy of all possible deexcitations in the molecule, and will be 
denoted by Cf^^^, (corr, deexc). Its form is slightly similar to the induction part: 



^'^'^S.'J^^^^'deexc) =Y,at 



{IaIa)0, 



{n-} 

X [(^s(A)|4j^B(^-}))®(^B({n-})|4j^B(A'))]^ (27) 

The summation in the above equation runs only over states of the molecule ipB{{n^}) with 
energy lower than the reference one, i.e. if C0n-,A = Ea — -E{n-} is positive, and hence Un-^A 
corresponds to the possible deexcitations of the molecule. This term does not have a simple 
physical interpretation, but as shown in Ref. [sil] it leads to a different QED retardation of 
the long-range potential than given by the classical Casimir-Polder formula [4lt]. Note also 
that without this extra term, the second-order interaction energy in the long-range could 
not be written correctly in terms of molecular properties of the isolated subsystems in the 
case when deexcitation may occur. Obviously, a similar term will be needed if an atom is 
in an exited state and molecule in its ground state. Then Cf^f^y (corr, deexc) would depend 
on the the dynamic polarizabilities of the molecule at frequency u corresponding to possible 
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atomic deexcitation. However, this holds only for atomic excited S states, as if the atomic 
state was P, D etc. then the whole formalism presented here would not be longer valid due 
to nonvanishing first-order energy in the multipole approximation. 



III. NUMERICAL RESULTS AND DISCUSSION 

We have applied the theory exposed above to the interactions of the ground state rubid- 
ium atom Rb(2S) with C0(3n), 0H(2n), NH(iA), and CH(2n), and of the ground state 
lithium atom Li(^S) with CH(^n). As discussed in the Introduction, these molecules have 
been successfully decelerated and are the best candidates for sympathetic cooling by colli- 
sions with ultracold rubidium atoms. At present no ab initio code allows for the calculations 
of all components of the dynamic polarizability tensor for open-shell linear molecules. There- 
fore, in our calculations we computed the polarizabilities appearing in Eqs. ( I23ll and (l22ll . 
from the sum-over-states expansion, Eq. ( IT6l) . The appropriate transitions moments to 
the excited states and excitation energies were calculated using linear response formalism 
with reference wave function obtained from the multireference selfconsistent field method 
(MCSCF) with large active spaces. For some electronic states the convergence of the sum 
in Eq. (fT6|) was not very fast, and we had to include over 100 states in the expansion. The 
DALTON program was used for linear response calculations. We have checked the conver- 
gence of the expansion (1161) by comparison of the static parallel components obtained from 
the MCSCF calculations with the results of finite-field restricted open-shell coupled cluster 
calculations with single, double, and noniterative triple excitations, RCCSD(T). The finite 



field calculations were done with the molpro code [43|]. Note parenthetically that finite field 
calculations can correctly reproduce the parallel component of the polarizability tensor, but 
fail for the perpendicular component due to the symmetry breaking. The nondiagonal com- 
ponents cannot be obtained from finite-field calculations. We have attempted to use the 
multiconfiguration interaction method restricted to single and double excitations (MRCI), 
but due to the convergence problems, we could get in this way only a few (up to ten) excited 
states. The Li and Rb atoms polarizabilities at imaginary frequency was taken from highly 



accurate relativistic calculations from the group of Derevianko |44|. 

In order to judge the quality of the computed long-range coefficients we have computed 
cuts through the potential energy surfaces of Rb-CII(^n) at a fixed distance -R=30 bohr 
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from the atom to the center of mass of the molecule. The zero of the angle 9 corresponds 
to the rubidium atom on the H side of CH. In these calculations we have employed the 
supermolecule method. The potential was computed as the difference, 



V 



TTiSM TTiSM TTiSM 

-^AB ~ -^A " -^B 5 



(28) 



where denotes the energy of the dimer computed using the supermolecule method SM, 
and E^, X=A or B, is the energy of the atom X. For the high-spin states (triplet for Rb(^S)- 
0H(2n), Rb(2S)-CH(2n), and Li(2S)-CH(2n), quartet for Rb(2S)-CO(3n), and doublet for 
Rb(2S)-NH(i A) we used the restricted open-shell coupled cluster method restricted to single, 
double, and noniterative triple excitations [RCCSD(T)]. Since the low-spin states have the 
same asymptotic behavior as the high-spin states there was no need to compute the ab 
initio poi nts explicitly. The RCCSD(T) calculations were done with the molpro suite of 
codes 431]. The distances in the diatomic molecules were fixed at their equilibrium values 
Te corresponding to the electronic state considered. For CO(^n) Ve was set equal to 2.279 
bohr, for OB.(^U) 1.834 bohr, for NH(iA) 1.954 bohr, and 2.116 bohr for CB.{^Il) The 
angle ^ = 0° corresponds to the linear geometris CH-Rb, CH-Li, OH-Rb, NH-Rb and 
CO-Rb. In order to mimic the scalar relativistic effects some electrons were described by 
pseudopotentials. For rubidium we took the ECP28MDF pseudopotential from the Stuttgart 



library j46|, and the spdfg quality basis set suggested in Ref. [46|. For the hght atoms 
(hydrogen, carbon, nitrogen, and oxygen) we used the aug-cc-pVQZ bases [47]. The full 
basis of the dimer was used in the supermolecule calculations and the Boys and Bernardi 



48 



scheme was used to correct for the basis-set superposition error 

Before going on with the discussion of the long-range interactions in the dimers, let us 
compare the diagonal static polarizabilities of CO(^n), OH(^n), NH(^A), and Cll(^n) with 
the literature data. In fact, the data are very scarce. For OH the most recent calculations of 
Spelsberg 2J| date back to 1999 (see also some older references 52h54||). For CH the only 
calculation we found in the literature is the 2007 paper by Manohar and Pal Q]. To our 
knowledge no data for the excited states of CO and NH were reported thus far. An ispection 
of Table [T] shows a relatively good agreement with the results of Spelsberg 24] for OH. The 
differences are of the order of a few percent, 5.5% for the parallel component, and 6.3% and 
3.3% for the perpendicular xx and yy components, respectively. For CH Manohar and Pal 



56| reported only the value of the parallel component obtained from the analytical second 
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derivative calculations with the Fock space multireference coupled cluster theory restricted 



56| is 



to single and double excitations. The agreement of our result with the value of Ref. 
remarkably good: the two results agree within 0.4%. 

The long-range coefficients C^'^ for the interactions of the ground state rubidium atom 
Rb(2S) with CO(=^n), 0H(2n), NH(iA), and CH(2n), and of the ground state lithium atom 
Li(^S) with CH(^n) are reported in Tables HTHVIi Also reported in these tables are the 
values of the induction and dispersion coefficients C^'^(ind) and C^'^(disp). First we note 
that for all systems and most of the coefficients the induction part is as important as the 
dispersion. This is not very surprising since the Li and Rb atoms are highly polarizable, 
and the molecules suitable for the Stark deceleration have large dipole moments. Note 
parenthetically that the coefficient C|^(ind) for the interactions of the 11 state molecules 
and (ind) for the interactions of the A state molecules vanish for symmetry reasons. The 
leading contribution to the anisotropy of the potentials in the long range, as measured by 
the ratio C|°/Cg°, is quite substantial for all the systems. The ratio C^^/Cg" ranges between 
0.26 for Li-CH and Rb-CH to 0.43 for Rb-OH. The difference in the anisotropy due to the 
presence of terms C^'^/R^ is relatively modest, since the coefficients are at least one 
order of magnitude smaller than C^°. 

Comparison of the long-range anisotropy of the potential energy surfaces for the singlet 
and triplet A' and A" states of Rb-CH(^n) computed from the mutlipole expansion up to 
and including R~^^ and from the supermolecule calculations is illustrated in Fig. [2l The 
agreement between the long-range and supermolecule results is relatively good, although 
small deviations can be observed. For the A' state the agreement is good at small angles, 
and slightly deteriorates for 6 around 100°. The same is true for the A" state, showing that 
our computed coefficients are not the perfect representation of the asymptotic expansion 
of the RCCSD(T) potential for this system. It should be stressed here that the long-range 
coefficients reported in the present paper do not describe the asymptotics of any potential 
obtained from supermolecule calculations, since for most of the supermolecule methods the 
long-range asymptotics is not known. See, e.g. Ref. [57| for a more detailed discussion of 
this point. However, the data reported in the present paper can be used in the fits of the 
potentials, or in the case of lack of ab initio points at large distances, to fix the long-range 
asymptotics with some switching function |58|. 

To illustrate the importance of the long-range coefficients with n > 6 in Fig. [2] we 
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report cuts through the potential energy surfaces of the least (Rb-CH) and most (Rb-OH) 
anisotropic systems for a fixed distance R = 30 bohr. An inspection of this figure shows 
that the contribution of the coefficients beyond ri = 6 is very important. For Rb-CH the 
terms qualitatively reproduce the anisotropy of the potential. This is not the case for 
Rb-OH. The inclusion of all terms up to n = 8 gives the correct picture of the anisotropy, 
and the and R~^^ contributions are of minor importance at this distance. It follows 
from the comparison of the RCCSD(T) results with the data computed from the asymptotic 
expansion, cf. Fig. |2l that the short-range exchange-repulsion effects are negligible at 
this distance. Thus, our illustration of Fig. [3] truUy demonstrates the importance of the 
R~^ and higher terms in the multipole expansion of the interaction energy. Obviously, the 
importance of the contributions beyond the Cq depends on the distance R, but our plot 
clearly shows that in the region of negligible exchange and overlap the contributions beyond 
Cq are important. 

No literature data are available for comparison, except for the long-range coefficients for 
Rb-OH obtained by Lara et al. 22| by fitting the CCSD(T) potential energy surfaces in 
the A' and A" symmetries at large distances to the functional form of Eq. ( |T8|) . The values 
of the long-range coefficients taken from Ref. [22i] are included in Table lllli The agreement 
between the two sets of the results is very reasonable. The isotropic coefficients agree 
within 7%. The discrepancies of the anisotropic coefficients are of the order of 10 to 15%. 
Such an agreement is satisfactory given the fact that the fitted values effectively account 
for the higher coefficients that could not be obtained from the fit. The only significant 
difference is in Cg^. Here the difference is as large as 37%, but this coefficient is small, and 
most probably could not be correctly reproduced from the fitting procedure. By contrast, 
the values of agree relatively well, within 10%. 



Hut son 



59| estimated the lowest dispersion coefficients Cq°, Cg°, and Cg^ for Rb-OH 



by using the best available data for the static polarizability of Rb and OH, and the Slater- 
Kirkwood rules. He obtained an isotropic Cq^ coefficient of 149.4 a.u., 30% off our value. 
Such an agreement is reasonable given all the simplifications of the Slater-Kirkwood ap- 
proach. The values of the anisotropic coefficients C|° and C|^, 4.3 and 0.9 a.u., respectively 
are four and three times smaller than those reported in the present paper, and thus unrealis- 
tic. This shows that the applicability of simple semiempirical rules to anisotropic interactions 
in open-shell complexes is of limited utility. We have performed a similar analysis for other 
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complexes considered in our paper, and came to similar conclusions. 



IV. SUMMARY AND CONCLUSIONS 

In the present paper we have formulated the theory of long-range interactions between a 
ground-state atom in an S state and a linear molecule in a degenerate state with a non-zero 
projection of the electronic orbital angular momentum. We have shown that the long- 
range coefficients describing the induction and dispersion interactions at large atom-diatom 
distances can be related to the first and second-order molecular properties. The final expres- 
sions for the long-range coefficients were written in terms of all components of the static and 
dynamic multipole polarizability tensor, including the nonadiagonal terms connecting states 
with the opposite projection of the electronic orbital angular momentum. It was also shown 
that for the interactions of molecules in excited states that are connected to the ground 
state by multipolar transition moments additional terms in the long-range induction energy 
appears. All these theoretical developments were illustrated with the numerical results for 
systems of interest for the sympathetic cooling experiments: interactions of the ground state 
Rb(2S) atom with 00(^0), 0H(2n), NH(iA), and CH(2n), and of the ground state Li(2S) 
atom with CH(^n). Our results for the static polarizabilities of the OH and CH molecules 



are in a good agreement with the ab initio data from other authors [2j, |56| . For all systems 
considered in the present paper the induction contribution to the long-range potential was 
found to be important. Also the anisotropy of the long-range interaction, as measured by the 
ratio C|°/Cg°, is substantial, while the anisotropy due to the is of modest importance. 
Relatively good agreement between the multipole-expanded and ab initio RCCSD(T) results 
was found. In the asymptotic region, where the exchange effects are negligible, terms 
with n > 8 are very important, and cannot be ne glec ted. For Rb-OH we could compare our 
results with the fit of ab initio RCCSD(T) points [22|. In general, relatively good agreement 
was found, except for the small Cg^ coefficient. It was also found that the Slater-Kirkwood 
rules for the anisotropic long-range coefficients fail in the case of open-shell monomers with 
spatial degeneracy. 
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TABLE I: Diagonal Cartesian components of the static dipole polarizabilities (in a.u.) for CO(^n), 
0H(2n), NH(iA), and CH(2n). 



co(3n) 0H(2n) NHi 


(^A) CH(2n) 


Reference 


^'^a,, 17.97 


8.29 11.23 


15.80 


present 




8.75 






Ref. [24] 








15.86 


Ref. [56] 


^'"a^^ 12.75 


5.99 


9.15 


14.32 


present 




6.37 






Ref. [24] 


' ^^yy 9.91 


7.31 


9.15 


11.81 


present 




7.55 






Ref. [24] 
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TABLE II: Long-range coefficients (in atomic units) for Rb— CO(^n). Cn' is the sum C'^(ind) + 
Cn '^(disp). The number in parentheses denotes the power of 10. 



L - 


> 


1 


2 


3 


4 


5 


6 


^6 


(ind) 1.187(2) 




1.187(2) 










6 


(disp) 3.797(2) 




6.400(1) 




























ind) 















6 


disp) 




-1.415(1) 










ck'^ 






-1.415(1) 












ind) 


4.646(2) 




3.103(2) 








7 


disp) 


1.045(3) 




-1.193(2) 








7 




1.470(3) 




1.913(2) 








^7 


(ind) 






4.027(1) 








7 


disp) 






-1.550(1) 








7 








2.477(1) 










find") 7 4(16^3'! 




6 165r3l 




1 423('3") 






8 


[disp) 3.144(4) 




6.753(3) 




2.401(2) 








3.884(4) 




1.292(4) 




1.663(3) 












3 168('2') 




1 136('2~l 






8 


disp) 




-4.200(2) 




-4.667(1) 






ck'^ 

8 






-1.021(2) 




6.695(1) 






^9 


ind) 


4.166(4) 




1 QQSM"^ 




7 102r3l 




9 


disp) 


1.226(5) 




-2.757(3) 




-2.947(3) 




ck° 




1.642(5) 




1.719(4) 




4.155(3) 






(ind) 






3.299(3) 




8.970(1) 






disp) 






-6.668(2) 




-2.376(1) 












2.632(3) 




6.594(1) 




^10 


(ind) 5.376(5) 




4.472(5) 




4.844(4) 




1.609(4) 


'-"10 


(disp) 9.360(5) 




7.679(5) 




4.014(4) 




-4.700(4) 




1.474(6) 




1.215(6) 




8.858(4) 




-3.091(4) 


riLi 


(ind) 




1.592(4) 




2.822(4) 




3.965(3) 


CLI 


disp) 




-4.121(4) 




-6.491(3) 




-1.094(3) 


CLI 
<-"10 






1.180(4) 




2.173(4) 




2.862(3) 
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TABLE III: Long-range coefficients (in atomic units) for Rb-OH(^n). Cn ' is the sum C^'"(ind) + 



-iL,K , 



Cn'^(disp). The number in parentheses denotes the power of 10. The values in square brackets 



are the results of Ref. 



22|. 



L - 


> 




1 


2 




3 


4 


5 


6 




(ind) 1.339(2) [1.33(2)] 
(disp) 2.154(2) [1.92(2)] 
3.494(2) [3.25(2)] 




1.339(2) 
1.654(1) 
1.505(2) 


[1.33(2)] 
[1.80(1)] 
[1.51(2)] 










r<L2 
r<L2 


(ind) 
disp) 








3.010(0) 
3.010(0) 


[1.90(0)] 
[1.90(0)] 










Cfo 


(ind) 
disp) 




9.460(2) 
2.679(2) 

1.214(3) [1.04(3)] 






6.306(2) 
1.058(2) 
7.365(2) [6.30(2)] 








Cf2 
Cf2 


(ind) 
disp) 










-4.040(1) 
4.800(0) 
-3.560(1) [-4.00(1)] 










(ind) 8.518(3) 
(disp) 1.626(4) 
2.487(4) 






8.188(3) 
2.772(3) 
1.096(4) 






2.117(3) 
0.274(0) 
2.117(3) 








(ind) 
disp) 






3.294(2) 
3.856(2) 
7.712(2) 






-4.442(1) 
1.397(1) 
-3.023(1) 








(ind) 
disp) 




7.136(4) 
2.865(4) 
1.000(5) 






4.188(4) 
1.211(4) 
5.399(4) 




5.980(3) 
6.591(2) 
6.641(3) 




nL2 


(ind) 
disp) 










-1.811(3) 
1.790(2) 

-1.632(3) 




1.273(2) 
3.633(1) 
1.637(2) 






(ind) 5.815(5) 
(disp) 3.085(5) 
8.900(5) 






5.738(5) 
2.046(5) 
7.784(5) 






1.457(5) 
2.241(4) 
1.681(5) 




8.853(3) 
1.705(3) 
1.056(4) 


nL2 

'-'10 
'^10 


(ind) 
disp) 






2.765(4) 
2.167(4) 
4.932(4) 






-1.320(2) 
3.751(2) 
2.431(2) 




6.332(2) 
8.513(1) 
7.183(2) 
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TABLE IV: Long-range coefficients (in atomic units) for Rb— NH(^A). Cn' is tfie sum C^'^(ind) + 
Cn'^(disp). Tfie number in parentlieses denotes the power of 10. 



L ^ 
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C|'"(ind) 1.120(2) 




1.120(2) 










C(f'"(disp) 2.849(2) 




2.094(1) 










Cr" :-!.l)()l)l2) 




i.:-!2!)(2) 










(ind) 


4.436(2) 




2.958(2) 








Cf°(disp) 


2.456(2) 




1.678(2) 
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4.636(2) 
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6.830(3) 




1.216(3) 






C|'°(disp) 2.370(4) 




2.611(3) 




3.100(2) 






2.981(4) 




9.441(3) 




1.526(3) 





















Ci'"'(disp) 








1.416(3) 






8 








1.416(3) 








3 ^24,(4) 




2.167(4) 




4 089(^3'! 




C^'O(disp) 


2.808(4) 




2.036(4) 




9.281(2) 




9 


6.132(4) 




4.204(4) 




5.017(3) 














-7.744(2) 




C|'4(disp) 










1.442(2) 














-6.302(2) 




Cf(0(md) 5.116(5) 




4.777(5) 




8.395(4) 




4.889(3) 


Cfo°(disp) 6.581(5) 




1.656(5) 




5.146(4) 




1.630(3) 


Cfo° 1.170(6) 




6.433(5) 




1.354(5) 




6.518(3) 


Cfo'(ind) 








4.640(4) 




1.128(3) 


Cf(f(disp) 








4.258(4) 




3.692(2) 


<^10 








8.898(4) 




1.497(3) 
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TABLE V: Long-range coefficients (in atomic units) for Rb-CH(^n). Cn is tfie sum C^'^(ind) + 
Cn'^(disp). Tfie number in parentlieses denotes the power of 10. 



L - 







1 


2 


3 


4 


5 


6 


^6 


(ind) 


9.654(1) 




9.654(1) 










f) 


(disp) 3.888(2) 




2.844(1) 










() 




1.8,-,:K2) 
















(ind) 

















6 


disp) 






-7.650(0) 










6 








-7.650(0) 










^7 


(ind) 




-2.974(2) 




-1.982(2) 








7 


disp) 




-3.874(1) 




3.254(2) 








7 






-3.361(2) 




1.272(2) 








^7 


(ind) 








2 ^fiQ^ 1 '1 

^ . c \ X ) 








7 


disp) 








-9.037(0) 








7 










1.665(1) 










(ind) 


6.731(3) 




3 341 r3l 




8 140('2") 






8 


(disp) 3.151(4) 




7.499(3) 




6.711(2) 






8 




3.824(4) 




1.084(4) 




1.485(3) 






^8 


(ind) 






1.170(2) 




-7 964(^0") 






8 


disp) 






-8.110(2) 




-5.038(1) 






8 








-6.940(2) 




-5.834(1) 








(ind) 




-2 130r41 




-1.349(4) 




3 089 r3l 




9 


disp) 




1.180(3) 




4.095(4) 




1.888(3) 




ck^ 






-2.012(4) 




2.746(4) 




4.977(3) 






(ind) 








1.108(3) 




-7.006(1) 






disp) 








-9.600(2) 




-2.901(1) 














1.012(3) 




-9.907(1) 






(ind) 


3.591(5) 




1.630(5) 




5.587(4) 




3.526(3) 




disp) 


8.679(5) 




8.055(5) 




8.953(4) 




8.714(3) 






1.227(6) 




9.685(5) 




1.454(5) 




1.224(4) 




(ind) 






1.170(3) 




8.307(2) 




7.917(2) 


riLi 


disp) 






1.238(3) 




3.143(2) 




-2.617(2) 


<^10 








2.408(3) 




1.145(3) 




5.300(2) 
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TABLE VI: Long-range coefficients (in atomic units) for Li— CH(^n). Cn is the sum C^'^(ind) + 
Cn'^(disp). The number in parentheses denotes the power of 10. 



L - 


> 


1 


2 


3 


4 


5 


6 


^6 


ind) 4.969(1) 




4.969(1) 










f) 


(disp) 2.034(2) 




1.521(1) 










() 


















find) 















6 


disp) 




-4.227(0) 










6 






-4.227(0) 










^7 


(ind) 


-1.531(2) 




-1.020(2) 










disp) 


-2.863(1) 




1.771(2) 












-1.817(2) 




7.502(2) 








^7 


ind) 






1.322(1) 








7 


disp) 






-5.003(0) 








7 








8.220(0) 










findl 1 Q44('3") 




5 030r2l 




6 6,34('2") 






8 


(disp) 1.125(4) 




3.670(3) 




3.547(2) 






8 


1.319(4) 




4.173(3) 




1.018(3) 






^8 


ind) 




6 n25Cl"l 




-4 loncn"! 






8 


disp) 




-3.344(2) 




-2.741(1) 






8 






-2.742(2) 




-3.151(1) 








ind) 


-4.721(3) 




-4 619('2"1 




1 590r,3l 




9 


disp) 


2.573(3) 




1.663(4) 




1.019(3) 




ck^ 




-2.147(3) 




1.619(4) 




2.609(3) 






(ind) 






1.604(2) 




-3.607(1) 






disp) 






-5.096(1) 




-1.568(1) 












1.095(2) 




-5.175(1) 






(ind) 6.133(4) 




1.983(4) 




1.337(4) 




1.815(3) 


•^10 


(disp) 9.074(5) 




3.187(5) 




3.728(4) 




4.704(3) 




9.687(5) 




3.385(5) 




5.065(4) 




6.519(3) 




(ind) 




1.424(3) 




5.150(2) 




5.621(2) 


riLi 


disp) 




3.582(2) 




2.407(2) 




-1.328(2) 


<^10 






1.782(3) 




7.557(2) 




4.293(2) 
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FIG. 1: Relative orientation of the molecule-fixed frame {xyz) to the space-fixed frame [XY Z) 
with its Z axis directed along the vector R joining the center of mass of the molecule BC and the 
atom A. The three atoms are in the plane axz = ctxz- 
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FIG. 2: Comparison of the long-range anisotropy of the potential energy surfaces for the doublet 
and quartet A' (solid line vs. circles) and A" (dotted line vs. crosses) states of Rb-CH(^n) 
computed from the mutlipole expansion up to and including R^^^ and from the supermolecule 
calculations. 



29 



Rb-CH R=30 bohr, A' 



Rb-CH R=30 bohr, A" 




60 80 100 120 140 160 180 
0[, 



Rb-OH R=30 bohr, A' 





FIG. 3: Importance of the R~^, R~^, and R~^^ contributions to the long-range anisotropy of the 
potential energy surfaces for the A' and A" states of Rb~CH(^n) and Rb-OH(^n) computed from 
the mutlipole expansion. 
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